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1 Introduction 

t-H ! 

J> ■ Factorization of linear partial differential operators (LPDOs) is a very well- 

0^ | studied problem and a lot of pure existence theorems are known. The only 

known constructive factorization algorithm - Beals-Kartashova (BK) factor- 
ization - is presented in pQ). Its comparison with Hensel descent which is 
sometimes regarded as constructive, is given in [2J, where the idea to use BK- 
| factorization for approximate factorization is also discussed. It originates in 

one of the most interesting features of BK-factorization: at the beginning all 
O j. the first-order factors are constructed and afterwards the factorization con- 



dition^) should be checked. This leads to the important application area - 
namely, numerical simulations which could be simplified substantially if in- 
stead of computation with one LPDE of order n we will be able to proceed 
computations with n LPDEs all of order 1. In numerical simulations it is 
^ ' not necessary to fulfill factorization conditions exactly but with some given 

^ ■ accuracy, which we call approximate factorization. 

The idea of the present paper is to look into the feasibility of solving 
problems of this kind using quantifier elimination by cylindrical algebraic 
decomposition [3J. In this paper we are going to apply this approach to a 
hyperbolic LPD0 of order 2 with polynomial coefficients. 

2 Hyperbolic LPDO of order 2 

A bivariate operator of second order has general form 

A 2 = a 2 odl + and x d y + a 02 dy + a w d x + a id y + a o (1) 
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and is factorizable [lj iff" 

f ujai + a i - £(2a 20 w + an) 1 wai + a i - £(2a 20 w + a n ) 

a o = £< r ■ S + ■ x 

[ 2cl2ouj + an J 2a 20 uj + an 

x Q2o(qqi - £(«20<^ + Q n)) + ( Q 20^ + aii)(aio - £020) ^ 
2a 2 ow + an 

Here coefficients ajj = a,ij(x,y) are functions on two variables x and y; 
w is a distinct root of the following polynomial V 2 {z) := a 2 o-z 2 + a\\z + 
O02> ^(w) = 0; and £ is a linear differential operator of the form C = 
d x - tody. 

Let us introduce a function of two variables x, y 

^ _ ^ f ^Qio + Qqi ~ £(2a 20 u; + an) | u;ai + a i - £(2a 20 a; + an) ^ 
\ 2a 2 o^ + an J 2a 2 ow + an 

Q2o(qqi - £(a20^ + an)) + (020^ + aii)(aio - £a 20 ) 
2a 20 w + an 

and rewrite factorization condition ([2]) as aoo = TZ. 

Now suppose that ([T]) is a hyperbolic operator in canonical form, i.e. 

H 2 = dl - dy + a 10 d x + aoi^ + a o (3) 

which corresponds to a 2 o(x,y) = 1, an(x,y) = 0, ao 2 (x,y) = —1. In this 
case we have two roots u± = 1 and u 2 = —1, and function 7Z takes a form 

- £ > ^ QQ1 1 + ( Q io ± Q oi) 2 



where "+" corresponds to u>\ = 1 and "— " corresponds to u> 2 = —1. We 
rewrite ([!]) is a slightly different form which will more convenient for further 
use: 

1Z = C{S} + S^ with 5= ((«io + «oi)Au, = l; (4) 

[(aio - a i)/2, u = -1. 

3 Polynomial coefficients 

Let us suppose that operator H 2 has polynomial coefficients dij and regard 
cases. 

3.1 Polynomials of first degree 

We have 3 polynomials a^- of first degree with two variables x, y: aoo(x, y) = 
6 3 ic + 6 2 y + 6i, a w (x,y) = c 3 x + c 2 y + ci, a i(ic, y) = d 3 a; + d 2 y + di, then 
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(1) For the first root u)\ = 1 we have C = d x — d y and 

aw + aoi = (c 3 + d 3 )x + (c 2 + d 2 )y + (ci + d\) = f 3 x + f 2 y + /i 

with /i = (ci+di), / 2 = (c2+d 2 ), /3 = (c 3 +d 3 ). Then £(ai +a i) = / 3 -/ 2 
and 



^(i) = h - h , C/gg + hv + /i) 2 
1 2 4 

(2) For the second root W2 = -1 we have C = d x + d y and 
«io - oqi = (C3 - d 3 )x + (c 2 - d 2 )y + (ci - di) = /i 3 z + h 2 y + fri 

with hi = (ci - di), /i 2 = (c 2 - d 2 ), /i 3 = (c 3 - d 3 ). Then £(a w - aoi) 
/13 — /i 2 and 

h 3 -h 2 (h 3 x + h 2 y + /ii) 2 



(5) 



ft 



(i) 



+ 



Remark 1. Notice that functions fti and 1Z 2 coincide symbolically: 

\2 



s 3 - s 2 (s 3 x + s 2 y + sx)' 



(6) 



(7) 



but of course, the form of s, as functions of coefficients Cj, bj, dj will be 
different. 

Remark 2. Factorization condition for the operator §3$) has now very simple 
form 



ft (1) = a 00 (x,y) 
which yields 



S3 - s 2 (s 3 x + s 2 y + si) 2 

o 1 : = 03^ + b 2 y + bi 



\ s 3 = 0, s 2 s 3 = 0, s 2 = 0, s 3 si = 26 3 , 
\s 2 si = 26 2 , a? + 2(s 3 - s 2 ) = 4&i 

For instance, in case of the second root this system of equations has form 

'(c 3 -d 3 ) 2 = 0, (ca - d 3 )(c 2 - d 2 ) = 0, (c 2 -d 2 ) 2 = 0, 
(c 3 - d 3 )(ci - di) = 26 3 , (c 2 - d 2 )(ci - di) = 26 2 , 
,( C1 - d x ) 2 + (2(c3 - d 3 ) - (c 2 - d 2 )) = 46a 

and its solution gives all exactly factorizable operators of this type: 

aoo(x,y) = (ci - di) 2 /4 
a 10 (x,y) = c 3 x + c 2 y + c\ 
_ aoi (x,y) = c 3 x + c 2 y + di 



3.2 Polynomials of second degree 

Now we have 3 polynomials a^- of second degree with two variables x, y: 

aoo(x, y) = b 6 x 2 + b 5 xy + 6 4 y 2 + b 3 x + b 2 y + h, 

aio(x, y) = cqx 2 + c 5 xy + c 4 y 2 + c 3 x + c 2 y + ci, 
a i(x,y) = d 6 x 2 + d 5 xy + d 4 y 2 + d 3 x + d 2 y + di, 

then 

1. For the first root wi = lwe have C = d x — d y and 

aio+aoi = (cQ+d 6 )x 2 +(c 5 +d 5 )xy+(c^+di)y 2 +(c3+d3)x+(c2+d 2 )y+(c 1 +d 1 ) 

= kx 2 + f 5 xy + / 4 y 2 + hx + f 2 y + fi 
with fi = (a + di), Mi = 1, 2, 6. Then 

£(ai + o i) = 2(/ 6 x - f A y) + f 5 (y - x) + / 3 - f 2 

and 

-,(2) _ 2(f 6 x - Uy) + h{y-x) + h- f 2 
*i - 2 + 

+ ^ (8) 

2. For the second root uo 2 = — 1 we have C = d x + d y and 

aio-ooi = (c 6 -d 6 )x 2 +(c5-d 5 )x|/+(c 4 -d 4 )y 2 +(c3-d3)x+(c2-d2)y+(ci-di) 

= /i 6 x 2 + h 5 xy + /i 4 y 2 + /i 3 x + h 2 y + hi 
with hi = (q — di), Vi = 1, 2, 6. Then 

£(aio - a i) = 2(/i 6 x - /i 4 y) + h 5 (y - x) + h 3 - h 2 

and 

^(2) _ 2Q 6 x - fo 4 g/) + fe 5 (y - x) + /i 3 - h 2 

(h 6 x 2 + h 5 xy + /i 4 j/ 2 + h 3 x + h 2 y + hi) 2 
+ 4 l y J 

As above we have in fact one function 

n (2) = 2 ( S 6 X - s 4 y) + s 5 (y - x) + s 3 - s 2 + 

| (s 6 x 2 + s 5 xy + s A y 2 + s 3 x + s 2 y + si) 2 



Remark 3. As direct corollaries of linear differentiation of polynomials one 
can conclude that for a polynomial any finite degree n, function TZ^ has 
form 

n— 1 

TZ {n) = T^ {k) + (other terms) 

k=l 

and is a polynomial of degree r with r < 2n, with necessary condition of 
exact factorization being deg(aoo) < deg(TZ^). 

4 Problem setting for QE 

We use standard formal language of elementary real algebra, that is, Tarski 
algebra [3] and formulate a first simple case of approximate factorization of 
LPDO as a quantifier elimination (QE) problem. Namely, we shall consider 
the case of a hyperbolic LPDO of order 2 with polynomial coefficients of the 
first order. In this case we have 

(1) three polynomials aij(x,y) of first degree in the two variables x,y; 
the coefficients of these polynomials are also variables: 

a 00 (x,y) = b 3 x + b 2 y + b 1 

aio(x,y) = c 3 x + c 2 y + ci 
a i(x,y) = d 3 x + d 2 y + di; 

(2) one function 

<r>(i)r ^ S 3-S2 . {s 3 x + s 2 y + si) 2 
K ( \x,y) = — - — + 

with Si given by (0) or by © , i.e. Sj = Cj+dj for the first root and S{ = Ci — di 
for the second root; 

(3) a constant e; 

(4) constants m and n, which define a bounded rectangular region in the 
plane: —m < x < m, —n < y < n. 

Remark 4. Notice that the special form of the factorization condi- 
tion allowed us to reduce the number of variables needed for this QE prob- 
lem. Initially we had 9 variables b%, b 2 , &i, c 3 , c 2 , ci, dg, d 2 , di, but in fact it 
is enough to consider only the 6 variables si, s 2 , S3, bi, b 2 , 63. 

With all this given, let us consider the quantified formula of elementary 
real algebra (jf = <f>*(bi, Sj) which asserts that "for all x and y in the bounded 
region — m < x < m, —n < y < n, we have — e < aoo(x, y) — 7Z^ (x, y) < e." 
We wish to eliminate the quantifiers from (p*(bi, Sj). More precisely, we wish 
to find a formula of elementary real algebra <p' = <f>'(bi, Sj), free of quantifiers, 
such that if (j>'(bi, Sj) is true then (p*(bi, Sj) is true. That is, we wish to find 
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conditions on the coefficients of the initial polynomials dij(x, y) which imply 
that the function TZ^(x, y) differs not too much from one these polynomials, 
namely aoo(x,y), throughout the bounded region — m < x < m, —n < y < 
n. 

5 Synopsis of QE by CAD 

Let A be a set of integral polynomials in x\,X2 ■ ■ ■ , x r , where r > 1. An 
A-invariant cylindrical algebraic decomposition (CAD) of R r , r-dimensional 
real space, is a decomposition D of R r into nonempty connected subsets 
called cells such that 

1. the cells of D are cylindrically arranged with respect to the variables 

2. every cell of D is a semialgebraic set (that is, a set defined by means 
of boolean combinations of polynomial equations and inequalities); and 

3. every polynomial in A is sign- invariant throughout each cell of D. 
The CAD algorithm as originally conceived (3J 0] has inputs and outputs 

as follows. Given such a set A of r-variate polynomials and a nonnegative 
integer / with f < r, the algorithm produces as its output a description 
of an A-invariant CAD D of R r , in which explicit semialgebraic defining 
formulas are provided only for the cells of the CAD Dj of R^ induced (that 
is, implicitly determined) by D. The description of D comprises lists of 
indices and sample points for the cells of D. (Every cell is assigned an index 
which indicates its position within the cylindrical structure of D.) 

The working of the original CAD algorithm can be summarized as fol- 
lows. If r = 1, an A-invariant CAD of R 1 is constructed directly, using 
polynomial real root isolation. If r > 1, then the algorithm computes a 
projection set P of (r — l)-variate polynomials (in ) such that 

any P- invariant CAD D' of R r 1 can be extended to a CAD D of R r . If 
/ = r we set /'«—/ — 1 and otherwise set /'<—/. Then the algorithm calls 
itself recursively on P and /' to get such a D' . Finally D' is extended to D. 
In order to produce semialgebraic defining formulas for the cells of Dt the 
algorithm must be used in a mode called augmented projection. 

Thus for r > 1, if we trace the algorithm, we see that it computes a 
first projection set P, eliminating x r , then computes the projection of P, 
eliminating x r -±, and so on, until the (r — l)-st projection set has been 
obtained, which is a set of polynomials in the variable x\ only. This is 
called the projection phase of the algorithm. The construction of a CAD of 
R 1 invariant with respect to the (r — l)-st projection set is called the base 
phase. The successive extensions of the CAD of R 1 to a CAD of R 2 , the 
CAD of R 2 to a CAD of R 3 , and so on, until an A-invariant cad of R r is 
obtained, constitute the extension phase of the algorithm. 

Now we consider the quantifier elimination ( QE) problem for the elemen- 
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tary theory of the reals: given a quantified formula (known as a QE problem 
instance) of elementary real algebra 

<P* = (Qf+lXf+l) ■ ■ ■ (QrX r )<f>(xi, ■ ■ ■ ,X r ) 

where <j) is a formula involving the variables x\, x^, ■ ■ ■ , x r which is free of 
quantifiers, find a formula <j)'(xi, . . . ,Xf), free of quantifiers, such that <fi' is 
equivalent to (f>* . The QE problem can be solved by constructing a certain 
CAD of R r . The method is described as follows. 

1. Extract from (p the list A of distinct non-zero r-variate polynomials 
occurring in (j). 

2. Construct lists S and I of sample points and cell indices, respectively, 
for an A-invariant CAD D of R r , together with a list F of semialgebraic 
defining formulas for the cells of the CAD Df of R/ induced by D. 

3. Using S, evaluate the truth value of (ft* in each cell of Df. (By 
construction of D, the truth value of (f>* is constant throughout each cell c 
of Df, hence can be determined by evaluating <f>* at the sample point of c.) 

4. Construct (/>'(xi, . . . , Xf) as the disjunction of the semialgebraic defin- 
ing formulas of those cells of Df for which the value of (j)* has been deter- 
mined to be true. 

The above algorithm solves any given particular instance of the QE prob- 
lem in principle. However the computing time of the algorithm grows steeply 
as the number r of variables occurring in the input formula <fi increases. 

Collins and Hong |7j introduced the method of partial CAD construction 
for QE. This method, named with the acronym QEPCAD, is based upon 
the simple observation that we can often solve a QE problem by means of a 
partially built CAD. The QEPCAD algorithm was originally implemented 
by Hong. A recent implementation, denoted by QEPCAD-B, contains im- 
provements by Brown, Collins, McCallum, and others - see 0- QEPCAD-B 
has solved a range of reasonably interesting problems for which the original 
QE algorithm takes too much time. Nevertheless the worst case computing 
time of QEPCAD-B remains large (that is, it depends doubly-exponentially 
on r). 

6 Application of QEPCAD to BK-factorization 

We consider only the first simple case of approximate factorization described 
in Section 4. Using the notation of Section 4, we suppose that e, m and n 
have been given specific constant values, say e = m = n = 1, and we consider 
the formula (f)*(bi,Sj) which asserts that 

(\te)(Vi/)[(M < 1 A|y| < 1) => \a 00 (x,y)-TZ^(x,y)\ < 1]. (11) 
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We wish to find a formula cp'(bi,Sj), free of quantifiers, such that (p'(bi,Sj) 
implies 4>*(bi, Sj). 

Remark 5. It would be of greatest interest to find the most general such 
4>'(bi,Sj) - that is, to find quantifier-free cj)'(bi,Sj) equivalent to 4>*(bi,Sj). 
But as we'll see it seems that the time and space resources needed to do this 
are prohibitive. We'll also see that it is not as time consuming, yet hopefully 
still of interest, to find quantifier-free conditions merely sufficient for <p* to 
be true. 

We attempted to find a solution to the above QE problem instance by 
running the program QEPCAD-B with the quantified formula [TT] (rewrit- 
ten so that the variables bi,Sj appear explicitly, and the denominator 4 is 
cleared from the right hand side of the implication, see [12] below) as its 
input. The variable ordering used was (S3, s%, si, 63, 62, &i, x,y). The com- 
puter used for this and subsequent experiments was a Sun server having 
a 292 MHz UltraSPARC rise processor. Forty megabytes of memory were 
made available for list processing. However the program ran out of memory 
after approximately one hour and forty minutes. The program was exe- 
cuting the projection phase of the algorithm when it stopped. The first 
three projection steps - that is, successive elimination of y, x and b\ - were 
complete. 

Increasing the amount of memory to eighty megabytes did not help - 
the program still ran out of memory during the fourth projection step (that 
is, during elimination of 62)- 

6.1 Searching for quantifier-free sufficient conditions 

Of course a very special, but completely trivial, quantifier-free sufficient 
condition for our QE problem instance is the formula 

<p'(bi, Sj) := [61 = A 6 2 = A 63 = A si = A s 2 = A s 3 = 0]. 

It could be of some interest to look for partial solutions to (that is, quantifier- 
free sufficient conditions for) our QE problem instance in which some but 
not all of the variables bi,Sj are equal to zero. For example, recall that 
the given quantified (jlip - after rewriting so that the variables bi, Sj appear 
explicitly and the denominator 4 is cleared from the right hand side of the 
implication - is: 

(Vz)(Vy)[(|x| < lA\y\ < 1) \ib 3 x+4b 2 y+4:b 1 -2(s 3 -S2)-(s3X+S2y+si) 2 \ < 4]. 

(12) 

Suppose that we put 62 = s 2 = in (fT2|) . We obtain: 

(\/x){Vy)[{\x\ < 1 A \y\ < 1) =► \4b 3 x + 4&i - 2s 3 - (s 3 x + si) 2 | < 4] 



8 



which is equivalent to: 

(Vx)[(|x| < 1) => |46 3 x + 46i - 2s 3 - (s 3 x + si) 2 | < 4], (13) 

which we shall denote by ip*(bi, Sj). 

The following theorem shows that a partial solution to the special QE 
problem instance ijj*(bi,Sj) (that is, a quantifier-free sufficient condition for 
ip*) leads to a partial solution to the QE problem instance <j>* (that is, a 
quantifier- free sufficient condition for <fr*). 

Theorem 1 Suppose that ip'(bi,Sj) is a quantifier-free formula, involving 
only bi,b 3 , s\, s 3 , which implies ip*(bi,Sj). Then the quantifier-free formula 
tp'(bi, Sj) A 62 = A S2 = implies 4>*(bi, Sj). 

► Let bi, Sj be real numbers. Assume tp'(bi, Sj) A 62 = A S2 = 0. Then 
i/j*(bi, Sj) A 62 = A S2 = is true, by hypothesis. Take real numbers x and 
y, with \x\ < 1 and \y\ < 1. Then 

\Ab 3 x+4b 2 y+Ab 1 -2(s 3 -s 2 )-(s 3 x+S2y+si) 2 \ = |46 3 x+46i-2s 3 -(s 3 a;+si) 2 | < 4, 
by virtue of (fT3|) (since \x\ < 1). Hence (fT2|) is true. ■ 

The above discussion suggests that it would be worthwhile to try to 
find a solution to the simplified, special QE problem instance ip* using the 
program QEPCAD-B. Putting (|13p into a slightly more general form, and 
hence reducing by 1 the number of variables in the formula, we obtain: 

(Vx)[(|x[ < 1) => \ax 2 + bx + c| < 4]. (14) 

A partial solution 9' (a, b, c) to (j!4[) could easily be transformed into a partial 
solution ^'(61,63,51,53) to ip* by setting a = —S3, 6 = 463 — 2S1S3 and 
c = 46i - 2s 3 - sj. 

We ran program QEPCAD-B with (|14p as its input. Eighty megabytes 
of memory were made available for list processing. After 191 seconds the 
program produced the following quantifier-free formula equivalent to (fT4"|) : 

c-b+a+4>=0/\c-b+a-4<=0/\ 
c+b+a+4>=0/\c+b+a-4<=0/\ 
[ 4 a c - b"2 + 16 a > \/ 4 a c - b~2 - 16 a > \/ 

[ b~2 - 16 a = A b~2 + 16 a > ] \/ [ b~2 - 16 a < /\ b - 2 a >= ] \/ 

[ b~2 - 16 a < A b + 2 a <= ] \/ [ b"2 - 16 a > A b + 2 a >= ] \/ 

[ b~2 - 16 a > A b - 2 a <= ] \/ 

[ b~2 - 16 a = A c - b + a + 4 > A c - b + a - 4 < ] ] . 

Since a = — s|, we have a < 0. We ran QEPCAD-B a second time, this time 
using the command 
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assume [a <= 0] . 

After 60 seconds the program produced the following somewhat simpler 
quantifier- free formula equivalent to (|14p under the assumption a < 0: 

c-b+a+4>=0/\c-b+a-4<=0/\ 

c+b+a+4>=0/\c+b+a-4<=0/\ 

[ 4 a c - b~2 - 16 a > \/ [ b > A b + 2 a >= ] \/ 

[ b < A b - 2 a <= ] \/ 

[ b~2 + 16 a = A c - b + a + 4 > A c - b + a - 4 < ] ] . 

It is possible to induce the program to produce an arguably even simpler 
solution formula using less computing time by making two separate runs of 
QEPCAD-B. The first run uses the command 

assume [a < 0] . 

After just 1.9 seconds the program produced the following quantifier- free 
formula equivalent to (|14|) under the assumption a < 0: 

c-b+a+4>=0/\c-b+a-4<=0/\ 
c+b+a+4>=0/\c+b+a-4<=0/\ 

[ b - 2 a <= \/ b + 2 a >= V 4 a c - b~2 - 16 a > ] . (15) 

The above formula is perhaps the most elegant and understandable of those 
obtained by applying QEPCAD-B to Formula For it is a slight improve- 
ment of (that is, slightly more compact than) a formula seen to be equivalent 
to it (under assumption a < 0) which is quite straightforward to derive by 
hand from (1140 using elementary properties of the parabola y = ax 2 + bx + c 
on the interval (—1,-1-1): 

[2a-b>=0/\a+b+c+4>=0/\a-b+c-4<=0]\/ 
[2 a + b >= A a - b + c + 4 >= A a + b + c -4 <= 0] \/ 
[2 a - b < /\ 2 a + b < /\ 4 a c - b"2 - 16 a > /\ 
a-b+c+4>=0/\a+b+c+4>=0]. (16) 

Remark 6. To derive by hand (16) from ()14[) under the assumption 
a < 0, one has to notice that function f(x) = ax 2 + bx + c has its max- 
imum value for f'(x) = 2ax + 6 = 0, that is, for x = —b/(2a), and con- 
sider three cases separately: (1) —b/(2a) < —1, (2) —b/{2a) > +1, and (3) 
— 1 < —b/(2a) < +1. For each of the above three cases one can then write 
down necessary and sufficient conditions for (j!4H to be true. For example, 
in Case 1, (fl4"|) is clearly equivalent to —4 <a + b + cf\a — 6 + c<4. After 
treating each of the above cases, we obtain (16) by forming the disjunction 
of the formulas corresponding to the cases. 

Of course, (15) for a < is not quite a complete solution to the QE 
problem instance of (|14p under assumption a < 0. To obtain a complete 
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solution we still needed to run QEPCAD a second time, this time for the 
case a = 0. For the second run we put a = in (|14p and use the command 

assume [b /= 0] . 

After 60 milliseconds the program produced the following formula equivalent 
to (|14p with a = under assumption 6^0: 

c -b+4>=0/\c-b-4<=0/\ 

c +b+4>=0/\c+b-4<=0 (17) 

This is immediately seen to be correct! Finally we could obtain a complete 
solution to (|14|) for a < by combining (15) for a < 0, (17) for b ^ a = and 
the formulae - 4<0/\c + 4>0 (for a = b = 0). In fact a simple and 
elegant way to achieve such a combination is to insert the disjunct a = 
into the last conjunct of (15): 

c-b+a+4>=0/\c-b+a-4<=0/\ 
c+b+a+4>=0/\c+b+a-4<=0/\ 
[ b - 2 a <= \/ b + 2 a >= \/ 

4 a c - b~2 - 16 a > \/ a = 0] . (18) 

7 Discussion 

As we remarked in Section 5 the worst case computing time of QEPCAD-B 
grows steeply as the number of variables in the given QE problem instance 
increases. Indeed, as is suggested by the results reported in Section 6, a 
complete solution of the QE problem instance (|12|) by QEPCAD-B using a 
reasonable amount of time and space seems to be unlikely for the foreseeable 
future. 

Nevertheless the results of Section 6 also suggest that QEPCAD-B could 
be of help in searching for certain kinds of sufficient conditions for (jlip . 
especially those which involve setting some of the variables to zero. 

We briefly mention here another kind of approach which a person could 
use to derive another kind of sufficient condition for (|11|) by hand. Namely, 
one could begin by expanding the polynomial aoo(x,y) — 1Z}~{x,y) in terms 
of x and y: 

a 00 (x,y) - TZ l {x,y) = (s 2 3 /4)x 2 + (-2s 3 s 2 /4)xy + (s 2 2 /4)y 2 + 
(b 3 - (s lS3 )/2)x + (b 2 - (s lS2 )/2)y + (h - (s 3 - s 2 )/2 - s\/A). 

By inspection of the terms on the right hand side of the above equation we 
see that a sufficient condition for (jlip is: 

I s_3~2 / 4 | < 1 / 6 A I 2 s_3 s_2 / 4 | < 1 / 6 A 

I s_2~2 / 4 | < 1 / 6 A lb_3 - (s_l s_3) / 2 | < 1 / 6 A 

I b_2 - (s_l s_2) / 2 | < 1 / 6 A 

I b_l - (s_3 - s_2) / 2 - s_l~2 / 4 | < 1 / 6. 
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The above sufficient condition is unlikely to be obtained in a reasonable 
amount of time and space using QEPCAD-B applied to (|12p . even if one 
issues assume commands. The number of variables involved is probably too 
big. However a version of QEPCAD-B which is planned for the future, which 
will have the capability to determine adjacency relationships amongst the 
cells of the partial CAD, could be of some use in analyzing certain topological 
properties of the truth set in six-dimensional space of the quantifier-free 
formula in bi,Sj above. 

Acknowledgements 

This paper has its origin in discussions between the authors at RISC-Linz 
during the second of S.M. to RISC-Linz in 2005. S.M. would like to thank 
Professors Franz Winkler and Bruno Buchberger, and all their colleagues 
and staff at RISC, for their hospitality during his stay. S.M. would also 
like to acknowledge helpful discussions and communications with Professors 
Daniel Lazard and Chris Brown. E.K. acknowledges the support of the 
Austrian Science Foundation (FWF) under projects SFB F013/F1304. 

References 

[1] R. Beals, E. Kartashova. "Constructively factoring linear partial differ- 
ential operators in two variables." TMPh 145(2): 1510-1523 (2005) 

[2] E. Kartashova, O. Rudenko. "Invariant Form of BK-factorization and 
its Applications." Proc. GIFT-2006: 225-241. Eds.: J. Calmet, W.M. 
Seiler, R.W. Tucker. Universitatsverlag Karlsruhe (2006) 

[3] G. E. Collins. "Quantifier elimination for real closed fields by cylindrical 
algebraic decomposition." LNCS 33: 134-183. Springer Verlag, Berlin 
(1975) 

[4] D. Arnon, G. Collins, S. McCallum. "Cylindrical algebraic decomposi- 
tion I: the basic algorithm." JSC, 13(4): 878-889, (1984) 

[5] C. Brown. "QEPCAD B: a program for computing with semi-algebraic 
sets using CADs." ACM SIGSAM Bulletin, 37(4): 97-108, (2003). 

[6] B. F. Caviness, J. R. Johnson (Eds.) Quantifier Elimination and 
Cylindrical Algebraic Decomposition. Springer Verlag, Berlin 
(1998) 

[7] G. E. Collins, H. Hong. "Partial cylindrical algebraic decomposition for 
quantifier elimination." JSC, 12(3): 299-328, (1991) 



12 



